TRANSIENT JETS IN THE SYMBIOTIC PROTOTYPE Z ANDROMEDAE 
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ABSTRACT 



OO 

■ We present the development of the collimated bipolar jets from the symbiotic pro- 

totype ZAnd that appeared and disappeared during its 2006 outburst. We monitored 
^ i the outburst with the optical high-resolution spectroscopy and multicolor UBVRq pho- 

^ \ tometry. In 2006 July ZAnd reached its historical maximum at C/ ~ 8.0. After ~lmag 

decline in mid- August, it kept its brightness at a high level of [/ ~ 9 up to 2007 Jan- 
uary. During this period, rapid photometric variations with Am ~ 0.06 mag on the 
timescale of hours developed. Simultaneously, high-velocity satellite components ap- 
peared on both sides of the Ha and H/3 emission line profiles. Their presence was 
transient, being detected to the end of 2006. They were launched asymmetrically with 
the red/blue velocity ratio of 1.2-1.3. From about mid-August onward they became 
symmetric at ~ it 1200 km s^^, reducing the velocity to ~ it 1100 km s^^ at their dis- 
appearance. Spectral properties of these satellite emissions indicated the ejection of 
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^2003-2006, visiting astronomer at the Department of Astronomy, Pennsylvania State University 



- 2 - 



bipolar jets collimated within an average opening angle of 6°.l. If the jets were expelled 
at the escape velocity then the mass of the accreting white dwarf is Mwd ~ 0.64 Mq. We 
estimated the average outflow rate via jets to Mjet ~ 2 x 10~^(i?jet/l AU)-^/^ yr~-^, 
during their August-September maximum, which corresponds to the emitting mass in 
jets, M?^ ~ 6 X 10-i°(i2jct/lAU)3/2 jv^0. During their lifetime, the jets released the 
total mass of M'^°^^'^ ~ 7.4 x 10^'^ Mq. Evolution in the rapid photometric variability 
and asymmetric ejection of jets around the optical maximum can be explained by a 
disruption of the inner parts of the disk caused by radiation-induced warping of the 
disk. 

Subject headings: stars: binaries: symbiotic - stars: individual: Z And - ISM: jets and 
outflows 



Introduction 



Collimated jets from astrophysical objects represent very exciting events. They have been 
detected from various kinds of objects. A common element in all jet-producing systems is the central 
accreting star, which implies t hat their accretion disks play an important role in the formation of 



this type of mass outflow (e.g. iLivio 



1997 



20041 . and references therein). 



Symbiotic stars are long-period interacting binaries (orbital periods are in the range of years) , 
in which a white dwarf (WD) accretes material from the wind of a cool giant. This process generates 
a very hot and luminous source of radiation (Th ~ 10^ K, Lh ~ 10^ — 10^ Lq) and accumulates 
a certain amount of material onto the WD's surface that enlarges its effective radius to about 
0.1 ii0 during their quiescent phases. During active phases a large, geometrically and optically 
thick circumstellar disk with radius of a few Rq develops around the accretor (jSkopal 1 12005| . 
Figs. 26 and 27). Therefore, the symbiotic stars are good candidates for producing high- velocity 
outflows in the form of collimated jets, mainly during outbursts. However , signatures of collimate d 



outflows were indicated for only 10 from roughly 200 known symbiotics (jBrocksopp 



et al. 



seven by resolving their spatial structure with the ra dio and X-ray imaging (e.g. I Taylor et al 



1986 



Crocker 



et al. 



Tomov et al. II2OO0I : ISchmid et al 



2001 



2004)1: 



2001 



Galloway k. Sokoloski Il2004l 'l and three by the optical spectroscopy (e.g. 



). Nevertheless, jets in symbiotic binaries are investigated 
very intensively from X-ray to radio wavelengths to understand better the p hysical process driving 



the collimated high-velocity mass outflow from a, c creting white dwarfs fe.g. ISokoloski et al. II2004I : 



Stute Camenzind II2005I : iKarovska et al. II2007I : iKellogg et al. II2007I : IStute Sahai Il2007l l 



Z And is considered a prototype symbiotic star. Its recent activity started from 2000 autumn 
and reached the optical maxima in 2000 December, 2004 September and 2006 July (Fig. 1). During 
the previous well observed large eruptions (1985, 2000), signatures of high- velocity outflows were 
detected through the b roadening of emission lines, P-Cygni type of profiles a.nd the extended wings 



of hydrogen lines (e.g. iFernandez-Castro et al. I Il995l : iTomov et al. 1 120031 . 12OO8I : ISokoloski et al 
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20061 : ISkopal et al. 1 120061 : ISkopal 1 120061 ) . but never in the form of collimated bipolar jets. The 
only signature for a nonsymmetric outflow from Z And was ind icated by the 5 GHz rad io map from 



2001 September, which showed a transient jet-like extension (IBrocksopp et al. 112004 1 



no counterpart in the optical spectrum could be r ecognized (IBurmeister fc: Leediarv 



However, 
Re- 



2003) 



cently, during the maximum of the 2006 outburst, ISkopal Pribulla I (|2006l ) discovered spectral 
signatures of bipolar jets from ZA nd. Their evidence in the optical spectra was confirmed by 



Burmeister fc Leedjarv I (j2007l ) and iTomov et al. I (120071 ). 



In this contribution we investigate evolution of the jet features from their first appearance 
in 2006 July to their disappearance at the end of 2006, determine their physical parameters, and 
discuss the disk-jet connection during the outburst. 



2. Observations and data reduction 

Our observations of ZAnd during its 2006 outburst were carried out at different observatories. 



(i) At the McDonald Observatory by the 9.2-m Hobby-Eberly telescope (jRamsev et al. Ill998l . 



HET in Table 1) just prior to the rise of the m ajor 2006 o utburst. Observations were performed 



using the echelle high dispersion spectrograph (jTull lllQQg ) with the 316g5936 cross-disperser, 3" 



fiber, and 2 CCDs (4096x4100) with 2x1 binning. Data were reduced using the standard IRAF pro- 
cedures, which involved bias, bad pixels, flatfield, scattered light, continuum, and cosmic rays cor- 
rections. Wavelength calibration was done using Th-Ar lamp and has an error of about 0.1 kms~^. 
Observations were recorded in 71 orders which cover the spectral regions 4060-5880, 6010-7760 A 
with the resolution of about R = 60 000 and S/N~200 in the continuum. For the purpose of 
this paper we use only regions containing Ha, H/?, He II 4686, He 1 4713 and the Raman-scattered 
O VI 6825 line. 

(ii) At the David Dunlap Observatory, University of Toronto (DDO) the high-resolution spec- 
troscopy was performed by the single dispersion slit spectrograph equipped with a Jobin Yovon 
Horiba CCD detector (2048x512 pixels of 13.5 /um size; thinned back illuminated chip) mounted 
at the Cassegrain focus of the 1.88-m telescope. The resolution power was 12 000 and 8 000 around 
the Ha and He II 4686— H/? region, respectively. Basic treatment of the spectra was done using 
the IRAF-package software. During each night two very different exposures were applied to obtain 
a well defined Ha profile and the continuum. Telluric sharp absorptions in the Ha region were 
eliminated with the aid of simultaneously observed standard star 50 Boo. 

(iii) At the Asiago Astrophysical Observatory (Asiago). Here the high-dispersion spectroscopy 
was secured by the REOSC echelle spectrograph equipped with a AIMO E2VCCD47-10 back illu- 
minated CCD detector (1100x1100 pixels of 13 /^m size) mounted at the Cassegrain focus of the 
1.82-m telescope at Mt. Ekar. The resolution power with the 200 //m width slit was approximately 
25 000 and 17000 at the Ha and He II 4686 - H/3 region, respectively. The spectroscopic data 
were processed with the MIDAS software package and a software developed at the Astronomical 
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Fig. 1. — Top panels show the U, B, V light curves (LCs) of Z And covering its two major eruptions that 
peaked in 2000 December and 2006 July. Arrows indicate times of spectroscopic observations. Absence of 
jets is denoted by "x". The lower block of panels displays spectral regions around Ha, H/3, and He II 4686 A 
during the 2006 outburst. The jet features are denoted by S~ and S+ and arrows. Dotted lines represent 
the continuum. 
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Table 1: Log of spectroscopic observations 



Date 


Julian Date 


Re 


gion 


T 

exp 


Obs. 


(dd/mm/yyyy) 


JD 24... 


(nm) 


(min) 




15/12/2005 


53720.127 


406 


- 776 


3x10 


HET 


04/06/2006 


53890.528 


462 


- 472 


30 


Asiago 




53890.551 


482 


- 492 


30 


Asiago 




53890.575 


652 


- 666 


30 


Asiago 


25/07/2006 


53942.809 


642 


- 671 


5 


DDO 


01/08/2006 


53949.837 


462 


- 492 


15 


DDO 


30/08/2006 


53978.688 


462 


- 492 


20 


DDO 




53978.706 


642 


- 671 


10 


DDO 


10/09/2006 


53989.602 


462 


- 472 


30 


Asiago 




53989.626 


482 


- 492 


30 


Asiago 




53989.650 


652 


- 666 


30 


Asiago 


18/10/2006 


54027.387 


641 


- 692 


33 


Ondfejov 


09/11/2006 


54049.418 


641 


- 692 


50 


Ondfejov 


23/11/2006 


54063.428 


641 


- 692 


50 


Ondfejov 


27/12/2006 


54096.905 


575 


- 690 


20 


OAO 




54096.919 


575 


- 690 


20 


OAO 


30/01/2007 


54131.250 


652 


- 666 


5 


Asiago 




54131.263 


652 


- 666 


30 


Asiago 




54131.286 


672 


- 686 


30 


Asiago 


11/04/2007 


54202.402 


642 


- 671 


5.5 


DDO 


22/06/2007 


54274.697 


642 


- 671 


6.0 


DDO 


31/07/2007 


54313.790 


642 


- 671 


3.0 


DDO 


27/08/2007 


54340.485 


652 


- 686 


30 


Asiago 



- 6 - 



Observatory of Capodimonte in Napoli. 

(iv) At the Ondfejov Observatory (Ondfejov) the high-resolution spectroscopy was performed 
by using the coude single-dispersion slit spectrograph of 2-m reflector and the BROR CCD camera 
with the SITe-005 800x2030 pixels chip. The resolution power at the Ha region was 10 000. Stan- 
dard initial reduction of CCD spectra (bias subtraction, flat-fielding and wavelength calibration) 
was carried out using modified MIDAS and IRAF packages. Final processing of the d ata was done 



with the aid of the SPEFO-package software developed at the Ondfejov Observatory (jHorn et al 



19961 : ISkodalll996l ) 



(v) At the Okayama Astrophysical Observatory (OAO) the high -dispersion spectra were se- 



cured with the high-dispersion echelle spectrograph (llzumiura I Il999l ) at the f/29 coude focus of 
the 1.88-m telescope. The dimension of the CCD (EEV 42-80) was 2048x4096 pixels of 13.5 /im^. 
The red cross-disperser was used. The spectral resolving power was 54 000 around the Ha region. 
The reduction and analysis was performed with the IRAF-package software. Two similar exposures 
were co-added to obtain a better signal-to-noise ratio (S/N). 

(vi) At the Skalnate Pleso and Stara Lesna (pavilion G2) Observatories, classical photoelectric 
UBVRq measurements wer e carried out by sing le-channel photometers mounted in the Cassegrain 



foci of 0.6-m reflectors fsee lSkopal et al. II2004I . in detail). The star BD+474192 (SAO 53150, V 
= 8.99, B - V = 0.41, U - B = 0.14, V - Rc = 0.10; CI in Fig. 3) was used as the standard 
star for both photoelectric and CCD observations. Fast CCD photometry was performed at Stara 
Lesna Observatory using the 0.5-m telescope (pavilion Gl). The SBIG STIO MXE CCD camera 
with the chip 2184x1472 pixels and the UBV{RI)c Johnson-Cousins filter set were mounted at 
the Newtonian focus. The size of the pixel is 6.8 yum and the scale 0.56"/pixel, corresponding to 
the field of view (FOV) of a CCD frame of about 24x16 arcmin. 

(vii) At the MonteBoo Observatory of the Masaryk University in Brno (Czech Republic), 
additional high-time-resolution optical photometry was obtained during two nights, 2006 Septem- 
ber 6 and 2006 September 11. The observational data were carried out using SBIG ST-8 dual 
chip CCD camera and Kron-Cousins BV{RI)c filter set attached at the Newtonian focus of the 
602/2780 telescope. The CCD chip has 1534x1020 pixels with the size of 9 //m and the scale 
2.04"/pixel (FOV~17.1xll.4arcmin). All CCD frames were dark-subtracted, flat-fielded and cor- 
rected for cosmic rays. Heliocentric correction was applied for all data points. CCD photometric 
data were processed with the MIDAS software package and a software C-MUNIPACK0, developed at 
the Masaryk University i n Brn o. Other details of the CCD photometric reduction were described 



by jParimucha &: Vafiko I (120051 ) . To verify the reality of Z And rapid variations during selected 
nights we measured additional comparison stars: C2 (GSC 03645-01592, V = 11.17, B — V = 0.53), 
C3 (GSC2.2 N012002338, B ~ 12.3 mag) and C4 (SAO 53133, V = 9.17, B-V = 1.36, U -B = 
1.11). The best results were achieved using the standard CI, because it is very similar to ZAnd 



^See the project homepage |http://c-munipack.sourceforge.net/| 
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Fig. 2. — Observed (full boxes) and modeled (solid thick li nes) SEDs at / around the 2006 maximum of 
Z And. Model for the red giant (dotted hne) was adopted from lSkopal I (j2005r ). Corresponding dates (1,2,3) 
are shown on the UBVRc LCs (•, A, +, □) inserted at the right-up corner. At times 2 and 3 the nebular 
contribution was negligible, while at the maximum (1) it increased by a factor of ^ 20 (denoted by N(l)). 
Radiation of a black body at the ionizing temperature of 90 000 K (bb) is scaled to a maximum of the hot 
object luminosity of ~ 10" (Sect. 3.3). 



during outbursts in both the brightness and the spectral type (see Figs. 2 and 3). The errors of 
the Z And - CI differences were mostly less than 0.01 mag. 

The journal of spectroscopic observations is given in Table 1. A correction for heliocentric 
velocity was applied to all spectra. Arbitrary flux units were converted to absolute fluxes with 
the aid of th e simultaneous UBV Rc photometry corrected for emission lines. The method was 
described by ISkopal et al. I (j2006l ). Uncertainties of such a continuum calibration are of a few 



percent for the star's brightness around 9 mag in the VR passbands (see ISkopal 1120071 . in detail). 
Observations w ^ere dereddened with -E'r-v = 0. 30 and resulting parameters were scaled to a distance 
of 1.5 kpc (e.g. iMikolajewska Kenyon Ill996l ). 



3. Analysis and results 
3.1. Photometric evolution around the maximum 



The top-right panel of Fig. 1 shows the UBV LCs of ZAnd covering its recent 2006 active 
phase. Prior to and after the peak July's magnitudes, the color index [/ — 1/ ~ (Figs. 1 and 2). In 
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July, during the maximum, the index had became temporarily much bluer {U — V ^ —0.55), which 
could be caused by a transient increase of the nebular emission. To understand better this behavior 
we modeled the U BVRcR .j flux points (F ig. 2 using the method of disentangling the composite 



continuum as described by ISkopal I (|2005l ). Our solutions correspond to a low temperature stellar 
source radiating at ~ 10 000 — 15 000 K that dominates the optical spectrum. This suggests a 
significant contribution from a warm shell. At the maximum (the model 1 in Fig. 2), the SED 
requires a larger nebular contribution (KitcUar /-^nebular = 2.3 in B, emission measure EM = 9.4 x 
10^°cm-3 at Te = 40 000 K), than prior to and after it (-Fsteiiar/i^nebuiar > 23 in B, EM ~ 1.3 - 
8.8 X 10^^ cm~^). The range of temperatures and scaling factors in the models of the stellar source 
corresponds to its effective radius = 12 ± ARq. We note that the absence of ultraviolet data 
did not allow us to determine the model parameters precisely. Uncertainties in and are as 
high as 30 and 50%, respectively. 

The effect of the transient increase of the nebular emission at the maximum could be interpreted 
as a result of a mass ejection from the active object into the particle bounded nebula, because new 
emitters will convert the excess of ionizing photons there into the nebular radiation. The effect of 
bri ghtening of symbiotic binary due to an increase in the mass- loss rate was originally suggested 



by iNussbaumer fc Vogel I (|l987l ) . 



In addition, we searched for a short-term photometric variability on the timescales of minutes 
to hours (Fig. 3). Just prior to and after the active stage we observe d an irregular variation 



within AB < 0.02 mag (top panels of Fig. 3; previously measured also by iGromadzki et al. 



Sokoloski Sz Bildsten Ill999l ). whereas around/after the maximum its amplitude increased by a factor 



2006 



of ~3, to AS ~ AV ~ 0.06 mag (bottom panels of Fig. 3; the first detection to date). Their source 
can be identified with the aid of the corresponding SED. The false (1 — 1.5) x 10^ K photosphere, 
which is associated with the disk, dominated the optical domain, while contributions from the giant 
and the nebula were negligible (Fig. 2). Therefore, the ~0.06mag variations had to be produced 
by the disk. Also the large size of the disk {Rd > R^ « 10 Rq) supports this interpretation ~ the 
timescale of the rapid photometric fluctuations is consistent with the dynamical time at the outer 
edge of the hot stellar source. 



3.2. Spectroscopic evolution 

Figures 1, 4 and 5 show evolution of the Ha, H/3, He II 4686, Hel4713A, and the Raman- 
scattered OVI 1032 profiles throughout the whole 2006-07 outburst. 

At the beginning of the 2006 active phase, prior to the initial rise (spectrum from 15/12/2005) 
as well as prior to the maximum (03/06/2006), the hydrogen and helium line profiles showed a simple 



^Rj fluxes wer e estimated from the measured _Rc magnitudes, corrected for the Ha emission, using transformations 
of Bessell I l|l983l ) for the less red stars. 
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Fig. 3. — Short-term variability in the optical continuum prior to the 2006-outburst (top left), after it (top 
right) and around its maximum (lower panels). Crosses represent individual CCD measurements, while full 
squares their means within 3.6 - 7.2 minutes. Comparison stars are denoted by CI, C2, C3 and C4 (see 
Sect. 2). 
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emission core. At these days, emission wings of Ha extended to ~ ±1500 and ~ itSOOkms"^, 
respectively, while those of H/3 persisted within ~ itSOOkms"^ on both spectra. The emission core 
of the He II 4686 line was placed symmetrically with respect to its reference wavelength (Fig. 1). Its 
strong flux in the 15/12/2005 spectrum, F(HeIl4686) = 2.3x 10~^^ erg cm~^ s~^, corresponded to 
the effective temperature of the ionizin g source, T^-'^- ~ 140 000 K (F(H/3) = 5.5 and F(Hel4471) = 



0.16x10 erg cm s : llijima Ill98ll ). In spite of a high temperature from nebular emission lines 



the Raman-scattered OVI 1032 line disappeared entirely (Fig. 5). 

Our first spectrum from the maximum (25/07/2006) revealed a complex of strong absorptions 
in the Ha profile. One was extended to ~ — lOOOkms"^ and others cut the emission core at 
~ — 220kms~^ and ~ — lOOkms"^ (Fig. 4). In addition, a pronounced S"*" satellite component 
was present and located at +1385kms~^. One week later (01/08/2006) strong absorptions at 
~ — 400kms~^ and ~ — 200kms~^ appeared in the H/3 and He 1 4713 A lines (Figs. 1 and 4). Two 
satellite components to the H/3 emission core were clearly recognizable and placed asymmetrically at 
~ — llOOkms^^ and ~ +1400 kms~^. From August 30th we observed typical broad emission wings 
with the extension of about it 2 000 and it 1 500 km s^^ in the Ha and H/9 profile, respectively. They 
were accompanied by pronounced bipolarly located satellite emissions. From October to December 
their fluxes were becoming fainter and at the end of 2007 January only some remnants in a form of 
faint spikes could be recognized (Fig. 4). We attribute these satellite components to bipolar jets. 
During the presence of jets, F(HeIl4686) fluxes decreased significantly (F(He II 4686) /F(H/3)^0.05) 



and their profiles were shifted blueward (Fig. 1) as during the 2000-01 maximurn (jSkopal et al. 



20061 ). The Raman line was not detectable (see also iBurmeister &: Leedjarv 1120071 ). It became to 



be visible on our 27/12/2006 spectrum as a faint redward-shifted emission around 6838 A (Fig. 5). 

Our spectra taken after the maximum, from 11/04/2007 to 27/08/2007, showed only a sim- 
ple emission core of Ha and H/3 lines with decreasing extension of their emission wings (Fig. 4). 
He II 4686 line became stronger, being located again around its reference wavelength as prior to the 
outburst (Fig. 1). The Raman emission band strengthened significantly and was comparable in the 
profile with that from quiescent phase (Fig. 5). 



3.3. A disk-like structure of the hot object and its evolution 



Observed SEDs from the optical maximum require a very low temperature of the hot object, 
Th ~ 1 — 1.5 X 10^ K (Sect. 3.1). In contrast, a significantly higher temperature of the ionizing source, 
T^^' ~ 140 000 and 90 000 K, was derived from t he hydrogen and helium emiss ion lines prior to the 



outbursts (Sect 3.2) and during its maximum ([Burmeister &: Leedjarv 1120071 ). respectively. Thus 
the hot active object was characterized by the two-temperature type of spectrum with <C T^'^' 
as during the 2000 maximum. This situation can be explained by a disk-like structure of the hot 
object viewed under a high inclination angle. Then the disk occults the hot ionizing source in 
the line of sight, w hile the nebula above/below the disk can easily be ionized (see Sect. 4.1. of 
Skopal et al. II2006I . in detail). 
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Heliocentric radial velocity [km s' ] 

Fig. 4. — Evolution in the Ha (H/3) broad wings along the 2006 outburst. The jet emission components 
(S^, S+, solid thin lines) developed during the optical maximum and were observable to the end of 2006 
December. We isolated them from the total profile by fitting its central emission with two Gaussian functions 

(dotted curves). The wind model of the broad Ha wings (Sect. 3.9) is shown here for clarity only to 
the spectrum from 15/12/2005 and 11/04/2007. Corresponding parameters are in Table 2. Fluxes are in 
10"-^^ ergcm"^ s~^ A~^. 



- 12 - 



Another effect resulting from the disk-structured hot object is the blueward shift of the 
He II 4686 A profile that develops during the outbursts, while out of active phas es it is placed 



symm etrically with respect to the reference wavelength (Fig. 1 here and Fig. 2 in ISkopal et al 



20061 ). The size of the He"*""*" zone was probably very small at the maximum (note the very low flux 
of the He II 4686 line during the maximum, Sect. 3.2) and thus the disk blocked a fraction of its 
redward shifted emission in the direction of the observer for the orbital inclination of about 75°. 



The disappearance of the Raman line during acti ve phases at high T^'^- (Sect. 3.2 and Fig. 5 
here; ISkopal et al. iBurmeister &: Leediarv 1120071 ) can also be a result of the disk-like structure 

of the hot object. The scattering process acts on the neutral atoms of hydrogen in the wind from 
the giant. So, if the ionizing source is indeed capable of produci ng the O"*"^ ions (note that t he 
ionization potential x(0"'"^) ~ 114 eV requires T^-'^- ~ 114 000 K (IMiirset fc Nussbaumer Ill994l )). 
bu t no Raman line is observed, implies that the O'*'^ zone is probably too small (cf. Appendix B 



m 



Skopal et alHbood ). so that the disk blocks its radiation in directions to the densest parts of 



the neutral wind at the orbital plane. However, the arise of a faint re(iw;ar(i-shifted Raman line 
on our spectra from 27/12/2006 and 30/01/2007 (Fig. 5) can be interpreted as due to the Raman 
scattering in the outer parts of the neutral giant's wind that move from the original O VI photons. 
This signals a gradual dilution of the disk. Following development of a strong Raman emission with 
a typical 6/Me-shifted shoulder in the profile (see Fig. 5) reflects total dilution of the optically thick 
material at the orbital plane, because the original O VI photons then can be scattered in the densest 
part of the neutral wind from the giant along the binary axis that moves against them (see also 



Schmid et al. Ill999l . for a general interpretation of the Raman line profile in symbiotic stars). The 



simultaneous disappearance of the disk and jets (cf . Figs. 4 and 5) is consistent with the necessity 
of the disk for the presence of jets (e.g. iLivio 1119971 ). The disk-jet connection is discussed below in 
Sect. 4. 

Finally, the disk- like structure of the hot active object in ZAnd was indicated for the first 
time by the model SEDs from the 1985 ou tburst, which revealed directly the presence of the two- 
temperature UV spectrum (Sect. 5.3.4. of ISkopal II2005I ). 



3.4. Measured parameters of the satellite components 

To isolate the jet components from the whole line profile, we formally fitted the emission line 
core and its extended foot component with two Gaussian functions. Then the residual jet emissions 
were compared with additional Gaussians. From their fitted parameters (the central wavelength, 
maximum, /, and the width a) we derived the jet's radial velocity RVs, flux Fa, = \p2j{ 1 a and 
the width FWUM<^ = 2y^21n(2)cr (Table 2). Corresponding fits are shown in Fig. 4. This 
approach and the resolution of our spectra allowed us to estimate uncertainties in the i?Vs values 
to 10-23 km s-i in Ha, FWRM^ widths were possible to adjust to observations within 0.2— 0.3 A 
and fluxes within 10—20% of the observed values. The smaller uncertainty corresponds to the 
stronger jet feature and vice versa. 
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Fig. 5. — Evolution of the Raman-scattered OVI 1032 line throughout the 2006 outburst. Top panel 
displays the U and V LCs with important periods in evolution of the Raman line, disk and jets. Bottom 
panel shows variation in the Raman line- profile along the outburst. Disappearance of the Raman line at high 
ionization temperature during active phases signals the presence of a large optically thick disk around the 
accretor (15/12/2005 - 23/11/2006). Its emergence (from 27/12/2006) indicates the dilution of the disk (see 
Sect. 3.3 in detail). Numbers at the left side represent shifts with respect to the level of the local continuum. 
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Our observations and those published by iBurmeister &: Leediarv I (|2007l ) and iTomov et al. 



( 20071 ) show an interesting behavior in RVs- From the jets launching in 2006 July to the beginning 
of 2006 August, the satellite components were located asymmetrically, being shifted by about 
+150kms~^ from the symmetric position (the top panel of Fig. 6). Their average velocity was 
1 280kms~^. Then, in the time of a few days, the jets became symmetric in their velocities. They 
had persisted at |i?Vs| ~ 1200kms^^ to about mid of September when reduced suddenly their 
\RVs\ to ~1 lOOkms^^ till the end of their detection. The middle panel of Fig. 6 plots the velocity 
ratio of the faster jet to its slower counterpart that reflects directly the evolution in the jets position. 
At the initial stage of the jets presence, this ratio was 1.2-1.3, whereas during the following period 
the velocity ratio was ~ 1. The bottom panel of Fig. 6 shows the flux ratio of the red jet component 
to the blue one, Fg+ /Fg- , for the Ha line as a function of the orbital phase. We interpret the 
decrease of this ratio, indicated around the inferior conjunction of the giant, as an occultation effect 
(Sect. 3.8.1). 

Relatively small width of the satellite components {FWHMs ~ 250kms~^, Table 2) with 
respect to, for example, the broad wings of the Ha line {FWZI{Ra) ~ 4 000kms-\ e.g. Fig. 4) 
implies that (i) the jets are highly collimated, and (ii) the emitting particles can be characterized 
with one constant velocity within the visible part of the jet. In the following two subsections we 
discuss some consequences of these jet properties. 



3.4-1- Collimated wind or jets? 



The winds of lumino us hot stars are accelerated by the radiation field, produced by the un- 
derlying photosphere (e.g. iLamers. &: Cassinelli Ill999l ). They are called line driven winds, because 
they are driven by absorption in spectral lines. In this process the photons transfer their energy into 
the kinetic energy of the wind, and, due to the Doppler shifts in the radial direction, the process 
acts throughout the entire wind, accelerating gradually its particles from the zero to the terminal 
velocity, foo- As a result, w e observe broad emission lines in the spectra of these stars. For symbi- 



otic binaries. ISkopal I (j2006l ) suggested that the broad wings of the Ha profiles {v^ 
are formed in the stellar wind from their hot components. 



> 2 000 km s" 



The relatively narrow width of the jet features suggests that they are driven by a different 
mechanisms. If the jets were driven like the wind, their observed profile would be a few thousands 
kms^^ broad, because of the large velocity dispersion in the radial direction (0 — Voo), and thus 
would hardly be detectable due to their low emission. Qualitatively, the narrow jet profiles imply 
that the jet engine works like a gun accelerating material through a short path (relatively to the 
jet length) to a high velocity, at which the shot-out-material continues to further distances. Also 
the connection between the disk and the jets, we discuss in Sect. 4, suggests that the jet launching 
mechanism is very different from that accelerating the wind. Generally, the jet engine transforms 
the kinetic energy of the inner parts of t h e disk to t he kinetic energy of th e flow in the vertical 
direction. In this respect. ISoker &: Regev I (j2003l ) and ISoker Lasota I (j2004l ) suggested a scenario 
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how jets could be blown from accretion disks. 



3.4-2. Geometry of the jet emitting region 

The narrow FWHMg allow us to assume a constant velocity of the emitting particles within 
the visible part of the jet. Then, if the broadening of the satellite components is caused by thermal 
motions, the corresponding Doppler width (~ 250kms~^, i.e. ~ 5.5 A in Ha) would require 
unrealistically high temperature of a few millions of Kelvins. In addition, the profiles of satellite 
lines have sharp edges, do not have typical wings, and their peaks are often flat (Fig. 1, e.g. 
10/09/2006). If the jet consisted of a parallel beam of particles, i.e. having the geometry of a 
narrow column, wc could not observe any real broadening. Therefore, wc ascribe the observed line 
width to the dispersion in the linc-of-sight velocity components of the jet particles. This can be 
satisfied if the jet emitting medium has the geometry of a narrow conus with the peak at the central 
object and characterized with a small opening angle ^o- 

Table 2: Parameters of Gaussian fits to the jet emission components: Radial velocity, -RVs [kms~^], flux, Fs 
[10~-^^ ergem"^ s~-^] and FWHMs [kms~-^]. The level of the local continuum, Fcont [10~-^^ ergcm"^ ^~^], 
opening angle 9o [°], emission measure of both jets, EMj^t [10^*cm~^] (Eq. 5), the mass loss rate through 
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3.5. Opening angle of the jets 



The conus geometry of the jet and the constant velocity of its particles (Sect. 3.4.2.) allow us to 
express the observed width of the satellite component as a function of its opening angle and the 
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Fig. 6. — Top: evoluti on of the jet radial velqcities, RVs - Large symbols rep resent our data (Table 2), 
while small ones are from^ Burmeister fc Leediarv ( 2007 ) and Tomov et al. ( 2007 ) for the Ha line. Measured 
values were shifted by ±1 lOOkms"^ to visualize better their relative position. Middle: RVs+ / RVg- ratios 
demonstrate the asymmetry in the jet radial velocities at the beginning of their launching. Compared is the 
ratio of opening angles (Sects. 3.5. and 4.2). Bottom: the ratio of the jet component fluxes, Fg+fFg-, for 
the Ha line as a function of the orbital phase. Their transient decrease around the phase 1.0 could be the 
effect of occultation of the red jet component by the giant (Sect. 3.8.1). 
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orbital inclination i. Assuming that the jets were launched with the velocity, Vjet, perpendicularly 
to the disk plane that coincides with the orbital one, the maximum dispersion of the line-of- 
sight jet velocity component corresponds to 2 x HWZIs = fjet cos(i — 6o/2) — Vjet cos(i + ^o/2) = 
2 sin(i) sin(0o/2), where the half width at the zero intensity of the jet (HWZIs) is in kms~^. Then 
for the measured central jet's velocity, RVs = v^^t cos(z), the opening angle can be approximated as 



2 sin 



HWZIs 
RVs tan(z) 



(1) 



A similar relation was also used by 
orbital inclination i = 76° (jSkopal 
as 



Shahbaz et al. I (|l997l ). Corresponding parameters from Table 2 



20oi) and adopting HWZIs = FWHMs, yield its average value 



6*0 = 6.°1±0.45 X Ai, 



(2) 



where the uncertainty represents a total differential of function (1) for uncertainties in FWHMs, 
RVs (Sect. 3.4) and that in the orbital inclination, Ai. 



3.6. The jet speed and the mass of the accretor 



It is generally considered that the jet velocity in all the jet-producing objects is of the order of 
the escape velocity from their cent ral stars, wh ich indicates that the jet-type outflo ws originate from 



the vi cinit y of the accretors (e.g. Livio 1 119971 . and references therein) . Recently, ISoker Regev 



(I2OO3I ) and lSoker fc Lasota I 

performed analytical estimates on how the plasma close to the 

central object, specifically in the boundary layer of the disk, could be accelerated to velocities larger 

than the local escape velocity. Our observations seem to be consistent with this scenario. The high 

velocities of jets, fjct ~ 5 000 kms~^, observed at large distances from the central star (Sec. 3.8.1), 

imply I'jet /''^escape > 1 for all possible distances of the jet ejection. However, from observations it is 

difficult to deter r nine a ccurately the distance from the accreting body, at which jets are launched. 

Soker &: Lasota I (j2004l ) also could not rule out the possibility that jets can be launched somewhat 

away from the accretor, up to ~ 3x the boundary layer radius. They noted that this case would 

require a strong perturbation in the disk. 

Therefore, for a rough estimate only, we will assume that the jets were launched from the 
WD surface at Vjet ~ ^^escape- Then the average value of Vjet = RVs/ cos(i) = 4960 ± 90 kms~^ 
{RVs = 1 200 ±22 kms"\ Table 2, i = 76°) and a typical WD radius of Rwd = 0.01 Rq correspond 
to the mass of the accretor in Z And, Mwd ~ 0.64 Mf;^. We note that this valu e satisfies perfectly 



the re lation between mass and radius of WDs (e.g. IShapiro &: Teukolsky Ill983l ). ISchmid &: Schild 
(119971 ) derived Mwd = 0.65 ± 0.28 Mq from the spectroscopic orbit for i = 47°. Our value of Mwd 
derived for i = 76° requires just a larger total system mass of ~ 3.2 Mq to satisfy the mass function 



of 0.024 Mq (jFekel et al. Il2000l ). However, in the real case, the same value of Mwd can be obtained 
for different combinations of the vjet/^escape (> 1) ratio and the distance of the jet ejection. 
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Finally, we note that the large effective radius of the active star during the optical maximum 
derived from the SED {R^ ~ 12 Rq ^> -RwDj Sect. 3.1) does not contradict to launching the jets 
from the vicinity of the WD surface, because of the disk-like structure of the hot active object seen 
under a high inclination angle (Sect. 3.3). 



3.7. Radiation of jets 



In this section we derive some constraints supporting that the jet emission is due to the 
photoionization of hydrogen and that the medium is optically thin in the direction of the observer. 



(i) According to Fig. 3 of iBurmeister &: Leediarv I (|2007l ) and the fluxes from Table 2, the 
average Balmer decrement -Ps (Ha) /Fg (H/3) ~ 2.4, which is close to the theoretical value of 2.75, 



given by the recombination process at the electron temperature of 20 000 K (e.g. iGurzadyan 1119971 ). 
Equivalent widths presented by Burmeister &: Leedjarv were converted to fluxes with the aid of our 
continuum calibration in Table 2. The somewhat lower values of the observed decrement than the 
theoretical one could be caused by a partial opacity of the nebular jet medium in the Balmer lines. 
Generally, the opacity in the Ha line is larger than in the H/3. The effect, however, should not be 
significant. Also the high orbital inclination and the small value of Oq correspond to a relatively 
small intersection, /, of the jet conus with the line of sight, and thus to a low value of the optical 
depth (r oc opacity x / ), which supports rather optically thin regime in the direction of the observer. 
For a comparison, the jets produced by the symbiotic star MWC 560 are optically thick, because 
they are seen nearly pole-on, i.e. their I and thus r are very larg e, which gives rise the blueshifted 



jet absorption components in the spectrum (ISchmid et al. Il200ll ) 



(ii) If the jet emission is due to the recombination of free electrons with protons, then it can 
be created only within the ionized fraction of the ejected material. In Appendix A we calculate 
the extension of the jet nebula as the distance from the central hot star, at which ionizing photons 
are completely consumed by jet particles along paths outward from the ionizing star. According 
to Eq. (A2) the jet radius is limited by its Stromgren sphere at the distance rg from the ionizing 
source, i.e. 



R 



jet = rs 



f 3Lph --2\ 



(3) 



where the parameters, Lph, OiBiH, Tg), n^^t are explained in Appendix A. In Section 3.8 and Figure 
7) we demonstrate that the jets radius, determined independently from their observed luminosity 
and the opening angle (Eq. (7)), agrees well with the radius of the Stromgren sphere, given by 
the luminosity and the temperature of the ionizing source. This result confirms that the satellite 
components are really created by recombinations within the jets, ionized by the central star. 

Evolution in the fluxes of jets is consistent with such an ionization structure. Between the 
beginning of 2006 A ugust and the end of September, their fluxes settle d (with some fluctuations) 



around a maximum (IBurmeister &: Leedjarv 



2007 



Tomov et al. 1120071 . this paper). No increasing 
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trend, as could be expected from an increase of the emission measure of jets (see Eq. (5)) due 
to their expansion, was observed during this period. This imphes that the emitting mass in jets 
was approximately constant, and thus suggests that the satellite components represent only the 
illuminated part of jets ionized by the hot central object, which is bounded by the Stromgren sphere. 
Figure 9 (see Appendix B) shows that the number of hydrogen ionizing photons is approximately 
constant for temperatures between about 50 000 and 130 000 K at a fixed luminosity of a black-body 
source. Therefore, the hot star produced roughly constant flux of ionizing photons, and thus also 
relevant fluxes of jets, in spi te that its temperature increa sed from 75 000 K to 115 000 K during 
the above mentioned period (jBurmeister &: Leediarv 1120071 ). Observed fluctuations were probably 
caused by those in the luminosity of the ionizing object. 



3.8. Mass loss through the jets 

The recombination process of the radiation by the jet plasma, the resulting fluxes, and the 
geometrical and kinematics parameters of the jets put some constraints to determine the mass loss 
rate through jets. Assuming that the jets were expelled into the solid angle Afi = 27r[l — cos{9q/2)] 
of the jet nozzle (Sect. 3.4.2.), then the corresponding mass-loss rate, Mjet, and the mean particle 
concentration, njet, in the jets volume, are connected via the mass continuity equation as 

Mjet = Afi Rf^^ njet-Ujet , (4) 

where /i is the mean molecular weight and mn is the mass of the hydrogen atom. Total luminosity 
produced by jets through the recombination transition of the Ha line, Ljet(Ha), is related to the 
line emissivity, e^nenp (ergcm~^ s~^), by 

Ljet(Ha) = Ea HeUpdV = nL Vjet = EaEM^et, (5) 

where Ea is the volume emission coefficient in Ha, Ue and np are concentrations of electrons and 
protons and Vjct and EMj^t are the volume and emission measure of the jets, respectively. For 
the optically thin medium of jets, the luminosity can be determined from the observed fluxes as 
Ljet = 4:iTcP X Fs- Further, Eq. (5) assumes a completely ionized medium (i.e. Ue = Up = njet), 
radiating at a constant electron temperature, i.e. is constant throughout the jet emitting volume. 
According to Sect. 3.4.2. we approximate its geometry by the conus, i.e. 

Vi,t = X (6) 

Substituting Eq. (6) into Eq. (5), we can express the jet radius by means of the parameters, obtained 
directly from observations (-Z>jet and AO), as 



3Ljet(Ha) 
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This expression allows us to rewrite Eq. (4) as 

[Any/' fL 



Mi 



jet 



jet \ ^^jct 



\njetj \[Lq] 

and/or function of the radius iijet as 



''jct -1 
— |- MqYT , 



(8) 



jet 



6 X 



( 



'jet ' -'^jet \ ^jet 



1/2 



[kms ^ 



(9) 



where the factor = 5.9 x 10 ^ or 1.1 x 10 ^ and ^2 = 3.6 x 10 ^ or 5.9 x 10 ^ for luminosities 
in Ha or H/3, respectivel y. Volume emissio n coefficients for Tc = 2 x 10'^ K, Eq, = 1.83 and = 
0.682 X lO-^^ergcm^s"^ dOsterbrock Ill989l '). 



3.8.1. Additional constraints for Mjet 

According to Eqs. (8) and (9), determination of Mjet requires to estimate reasonable values 
of njet or i?jet- Figure 7 plots dependencies between these parameters. Below we discuss some 
constraints that allow us to determine more accurate ranges of the jet parameters. 



(i) ISoker Lasota I (|2004l ) introduced conditions for thermally launching jets from accretion 
disks around WDs. In their model the accreted material is strongly shocked due to large gradients 
of physical quantities in the boundary layer and cools on the timescale longer than its ejection 
time from the disk. The model requires large accretion rates of Mace ^ lO^^MQyr^^. Some 
observational evidences of the disk-jet connection, as described in Section 4, could support this 
scenario. In agreement with this suggestion and the fact that there is a large disk encompassing 
the WD (Sects. 3.1. and 3.3.), we can exclude rates of Mjet < 5.5 x lO^'^Moyr"! (i.e. Mace < 
5.5 X 10-6 Mgyr-^ Sect. 4.1) that correspond to njet > 10^^ cm-^ and Rjet < 0.086 AU (see 
Fig. 7), because of too small jet radii that would be significantly occulted by the disk for the orbital 
inclination of ~ 75°. Observations do not indicate this case (e.g. Table 2, Fig. 6). From the other 
side, we exclude rates Mjet > 5.5 x 10"^ M© yr"! (i.e. njet < 10^ cm i?jet ^ 9AU, Fig. 7), because 
this would require too high accretion rates of > 5.5 x 10""^ Mq yr^^ (see also Sect. 4.1 ). As the mass 
of the disk can be estimated to ~ lO""^ - 10"^ Mq (e.g. Rn = 20 Rq, H/Rb = 0.3 (jskooal IboO^ ) 
and average density log p = -7 to -8), the accretion rates of about 10~^ Mq yr~^ would exhaust the 
disk during a short time. However, observations indicate the presence of a large disk with jets for 
about six months, from 2006 July to 2007 January (Sect. 3.3, Fig. 5). These conditions thus limit 
Mjet to 5.5 X 10~6 > Mjet > 5.5 x lO"''M0yr"^ and consequently, 10^ < njet < 10^^ cm~^ and 
0.1 < iijet < 9AU. 

(ii) During active phases of ZAnd the emission measu re produced b y the ionized hot star 



wind was determined to EMy^ ~ 3 x 10 cm (Table 1 in ISkopal II2006I ). It represents the so- 



called low-temperature nebula, which is subject to eclipses in active symbiotic systems (ISkopal 
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20051 . |2006| ). and thus it is located within the radius of the giant, Rq, around the hot component. 
For ZAnd, the average particle concentration in the wind around the hot star within the radius 
Rq = 106 Rq is nw = {EM^^ /VgY^'^ = 1-3 x 10^*^ cm~^. For a comparison, the mean concentration 
of the jet particles within this distance, njet(0.5AU) ~ 7.5 x lO^'^cm"'^, which corresponds to 
Mjet ~ 1.3 X lO~^M0yr~^ (Fig- 7). This suggests that n-^et > "-w and also the normalized mass 
loss rate (e.g. into 1 sr) via the jets is larger than that through the wind. 



(iii) The transient decrease in the flux ratio, -Fg+/-Fg- < 1, measured during the inferior con- 
junction of the giant (Fig. 6), could be ca used by the occultat ion of the red jet component by 
the stellar disk of the giant (also noted bv lSkopal &: Wolf 1120061 ). As 2Rq ~ 1 AU, this observa- 
tion suggests that the radii of jets are <; 1 AU. However, uncertainties i n the orbital elements do 



not allow a niore accurate estirnate (s ee iFormiggini &: Leibowitz 1 1 19941 : iMikolajewska Kenyon 



1996 



Skopal I Il998l : iFekel et al. 1 12000) . Finally, we note that some theoreti cal works, devoted to 



modeling the jets in symbi otic stars, also adopt the jet radius of 1 AU (e.g. IStute Camenzind 



2005 



Stute &: Sahai 1120071 ). Therefore, we determined the mass loss rate Mjet in Table 2 for a 
representative value of iijct = 1 AU. 



3.8.2. Emitting mass in jets 
The mean concentration and the volume of jets determine their emitting mass as 

MjX =/^mHnjetl^et, (10) 

which can be rewritten with the aid of Eqs. (6) and (7) as 



where ^3 = 5.7 x 10^^ or 9.3x 10~^ for luminosities in Ha or H/3, respectively. The average value of 
9q = 6.°1 (Sect. 3.5), Ljet = 0.32 Lq in Ha (from the 2006 August ~ September maxi mum) and R-^pt 
1 AU yield the emitting mass of = 6 x lO^^*' Mq in both jet components. iTomov et al. 



( 20071 ) derived significantly larger value of = 1.4 x 10 Mq for the similar luminosity of the 



jet 

jet pair they observed on 12/08/2006. This difference results from their adopted mean density of 
only lO^cm"^, which requires very large size of one jet, i?jet=40AU. However, such an extension 
of the jet is beyond the limits we estimated in Sect. 3.8. 



3.8.3. Total mass released hy jets 



With respect to our approximation, in which the jet medium has a constant particle density 
and emits within its Stromgren radius, we approximate the total mass released by the jets during 
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the time of their detection, ATjet, as 

Mj*°*"' « Mjet X ATjet, (12) 

where Mjet is the average mass loss rate via jets. Values from Table 2 (Mjet = 1-75 x 10^^ Mq yr~^ 
for the jet pair and ATjet = 154 days) yield M*gf ^' 7.4 x 10"^ Mq. 



3.9. Mass loss through the vi^ind 

It is generally well known that during active phases the hot components in symbiotic binaries 
lose their mass in the form of wind. To estimate the rate of the mas s loss produced by the active 



star through the stellar wind, we used the Ha method as suggested bv lSkopal I (|2006l ). The method 
assumes that the broad wings are due to kinematics of the photoionized and optically thin stellar 
wind. Profiles from the maximum were modeled with the same input parameters as those from the 
2000-maximum. The disk radius, i?D = 6.4 i?©, and the disk thickness at its edge, H = 1.92 Rq, 
were adjusted to the effective radius of the hot object derived from the SED assuming a flared 
disk with H/Rd = 0.3. Prior to and after the outburst we used i?D = 3.5 i?0 and H = 1.05 Rq. 
The results do not depend critically on these parameters. Examples of comparison between the 
modeled and observed profiles are shown in Figs. 4 and 8. Synthetic profiles fit well the observed 
wings for \RV\ > 300 km s^^. Corresponding mass loss rates through the wind, Mw are around 
2 X 10~^ Mq yr^^ (Table 2). Other model parameters are the terminal velocity, whic h runs betwee n 
1 500 kms^^ and 3000 kms^^, and the acceleration parameter /3 ~ 1.7 (cf. Eq. (2) in lSkopal II2006I ). 



It is of interest to note that the observation made just prior to the outbursts (the HET 
spectrum) revealed a more complex structure in the wind outflow. The O — C plot shows a faint 
and very broad (FWHM « 1 200 km s"^) emission component superposed on the extended blue Ha 
wing (Fig. 8). 



4. A disk-jet connection 



4.1. Observational evidence 



The red giant in ZAnd loses mass via the wind at Mq ~ 7 x 10 ^ Mq yr ^ ( Skopal 



200, 



a 



Table 3). According to the recent study on the wind accretion in binary stars (jNagae et al. II2004I ) 
the WD companion can accrete with efficiency up to 10%, i.e. at M™™'^ ~ 7 x 10~^ Mq yr~^ during 
quiescent phases of ZAnd (the total system mass of 3.2 Mq and the wind velocity of 30kms~^ 
at the accretor were adopted). This quantity corresponds to the accretion luminosity of ~ 70 Lq 
(parameters from Sect. 3.6), which is at le ast a factor of 30 smaller than the observed luminosity 
of the hot stellar object (e.g. ISkopal II2005I ). This discrepancy led to suggestion that an additional 
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Fig. 7. — Dependencies of i?jct (Eq. (7)) and Mjct (Eq. (8)) on njct- They were calculated for average values 
of Ljet = 0.32 Lq (to the end of 2006 Sept.) and 6^0 = 6.°1 (i.e. An = 8.9 x lO-^sr). Mjot here represents 
the mass loss rate by both the jets. C ompared are radi i of th e Stromgren sphere rs (Eq. (3)) for the hot 
star luminosity of 10'' L©, estimated by Sokoloski et al. ( 20061 ) (lower line), and for 2.2 x 10^ L© as results 
for scaling this quantity to our parameters of d — 1.5 kpc and £'b-v — 0.30 (upper line). A good agreement 
between both i?jet and rs radii confirms the recombination process as responsible for the radiation of the jet 
plasma, ionized by the central star. 




-2000 -1000 1000 2000 -2000 -1000 1000 2000 

Heliocentric radial velocity [km s"'] 

Fig. 8. — Two examples of the observed and synthetic Ha profiles from the ionized stellar wind of the hot 
star. The O — C plot demonstrates a good fit for |i?V^| > 300 kms^^. The models correspond to the mass 
loss rate of 1.3 and 2.5x lO-^Moyr-i for 15/12/2005 and 30/08/2006, respectively, Voo = 2 700kms-i and 
P = 1.7 (Sect. 3.9). Fluxes are in lO-^^ gj-gcm-^ s'^ A^^. 
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source of th e energy in a form of steadily b urning material on the white dwarf surface has to be 
present (e.g. iMikolaiewska &: Kenyon Ill992l ). 



During the recent active phase, the Mjet rate put limits for the accretion rate through the 
disk, Mace, be cause both the rate s are proportio i ial. Typical values of Mpt/Mam are in th e range 
0.01-0.3 fe.g. iLizano et al. lllQSsI : IPrinde \\l99^ : \Uvio \\199^ ). although llJvio et al. I (|2003l ) found 



that Mjet can almost be equal to M^cc , if the accretion energy is converted efficiently into magnetic 
energy and is emitted in the form of a magnetically dominated outflow or jet. Here we will consider 
the ratio Mjet/Macc ~ 0.1. Accordingly, the observed Mjet ~ 10"6(i?jct/l AU)^/^ jv/q yr-i thus 
requires M^cc ~ 10~^MQyr~^. Two possible mechanisms how to increase the accretion rate are 
either a disk instability or an increase in the mass loss rate from the giant. If the latter possibility 
had been the case, the normal red giant in Z And would have increased its Mq by 3 orders of 
magnitude - from a few times 10"'^ Mq yr~^ during quiescence to « 10~^ Mq yr~^ for a short time 
during the activity, to power the observed jets. This seems to be unlikely. Thus, assuming that Mq 
is more or less constant, a disk instability had to be responsible for the transient, but significant, 
increase in Mace through the disk. We summarize critical observations supporting this suggestion 
as follows: 



1. The rapid photometric variability originates in the disk (Sects. 3.1 and 3.3). Figure 3 shows 
the change from an irregular low-amplitude variation prior to the outburst to a slower, but 
higher-magni tude variation, obser v ed du ring the maximum when the jets were launched. 
According to ISokoloski &: Kenyon I (|2003l ). who examined the disk-jet connection in CHCyg, 
the disappearance of the fastest variations indicates that the innermost disk was disrupted 
around the time when the jets were produced. This interpretation is based on the relationship 
between the disk radius and variabilit y timescale (slower variatio ns come from larger disk 
radii and vice versa, as summarized by lSokoloski fc Kenyon II2003I ). Also the gas-dynamical 
modeling of the flow structures during the Z And outburst suggest ed that the sudden in crease 



of the accretion rate could result from the disruption of the disk (iBisikalo et al. 1120061 ). The 



disk disruption, i.e. the removal of material from its inner parts, implies that the jets were 
accretion-powered. 

Dramatic spectroscopic evolution around the optical maximum (2006 July/ August), as in- 
dicated by complex absorption/emission profiles accompanied by the asymmetric ejection of 
jets (Fig. 4, Sect. 3.2 here and Fig. 2 in iTomov et al. 1120071 ). could also be a direct result 
of the disk disruption. It is of interest to note that similar behavior was also observed dur- 
ing CHCyg outburst s in 1982-84, 1992-9 4 and 1998-00, but on somewhat longer timescales 
(see Figs. 1 and 4 of lskoual et'allbood'l. All the outbursts were followed by jet-extended 
features on the radio images ( Crocker et al. II2OO1I ). However, for ZAnd these features have 
been record ed for the first time. For example, during it s 2000-03 major outburst no jets were 



present (see ISokoloski et al. 



2006 



Skopal et al. II2OO6I ). 



3. Simultaneous presence of both the disk and jets during the period 2006 July - December and 
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their disappearance from 2007 January (Sect. 3.3, Figs. 4 and 5) suggest that bulk of the 
accretion energy was released in the form of jets. As a result the hot stellar source in Z And 
had reduced its radiation significantly and therefore we observed a rapid decline in the star's 
brightness from 2007 January (see Figs. 1 and 5). From this point of view, the absence of 
jets during the previous (2000-03) major outburst is consistent with a slower decline in the 
bright ness from the 2000- maximum (see Fig. 1) and the disk presence for a longer time (~1.5 



Skopal et al. II2006I ). In both the outbursts the presence of a large disk was terminated 



years, 

observationally by a small optica l rebrightening acco mpanied by the emergence of the Raman 



line (Fig. 5 here and Sect. 3.6 in lSkopal et al. II2006I ) 



We note that the disk-jet connection in Z And, as summarized in the points ab ove, has some 



simila rities with a hard state of the microquasa r GRS 1915+105 as described by iBelloni et al 



(|2000l ) and later di scussed by 



Livio et al. 



was pointed out bv lSokoloski &: Kenvon I (|2003l ) 



(120031 ) . For the symbiotic star CH Cyg this similarity 



4.2. Variation of jets and disk — a self-induced warping of the disk? 



The asymmetric jet episode in Z And lasted for less than one month - from the jets launching 
in 2006 July to the beginning of 2006 August (Fig. 6, Sect. 3.4). The velocity asymmetry in jets 
was so far revealed only for other objects and at significantly larger ti r nescales of years to decades 



(e.g. iHirth et al. Ill994l : IWoitas et al. 



2002 



Lopez-Martin et al. ll2003l : lNamouni 1120071 ). Thus, the 



very short duration of the asymmetric jets in Z And was probably connected with a disk instability 
induced during the optical maximum. Qualitatively, the instability could cause different conditions 
on the accretor's poles for the ejection of jets. Particularly, the anti-correlation of the opening 
angles to the jet velocities at the initial stage of the jets launching (Fig. 6) could be a possible 
cause of the observed jet's asymmetry. The larger the jet nozzle is, the lower ejection velocity 
can be expe cted to drive the j et flow at the same rate (Eq. (4)). This possibility was already 
suggested bv lHirth et al. I (|1994| ) for magnetically driven jets if the opening angles are different due 
to different pressure gradients. The short-term. Am ~ 0.06 mag, photometric variations represent 
another type of variability, which developed during the jet ejection (Fig. 3, Sect. 3.1). Their source 
was associated with the disk (Sect. 3.1). In the previous section we ascribed the evolution in the 
short-term variability to disruption of the inner parts of the disk. 

The origin of both these effects could become better understood by investigating a disk insta- 
bil ity due to an i ncrease of the hot component luminosity at the outburst maximum. According 



to iPringle I (|l996l ) irradiation of the disk by the central star can lead an originally planar disk to 
become twisted and tilted out of the orbital plane. In our case, additional energy liberated during 
the outburst by thermonuclear burning on the WD s urface increases formally the efficiency e of 
the accretion process (see Eq. (3.10) in |Pringle1ll997l ). which then could lead to radiation-driven 
warping close to the disk's center even for WD accretors. For Z And, the radiation-induced warping 
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occurs at all radii R > 0.6 Rq for the luminosity of ~ 10^ Lq around the optical maximum (e.g. 
Sokoloski et"al~ll2006l ^. the mass accretion rate from the wind of ~ 7 x 10 ^ M^yv ^ (Sect. 4.1) 
and the WD mass of 0.64 (Sect. 3.6). With respect to the disk radius of >IQRq (Sect. 3.1), 
the warping will act from inner parts of the disk. 

However, connections between jets and disk variability in symbiotics have not been studied 
theoretically in detail yet. Current theories on the warping and wobbling disks were elaborated 
for different types of objects and are ch aracterized by significantly larger tirnescales than those 
sugge sted by observations of Z And (e.g. iPringle I Il996l : iLivio &: Pringle 1 119971 : IWijers &: Pringle 



19991 ) 



5. Summary 

The main results of this paper can be summarized as follows: 

(i) Between 2006 July 19 and 28 our photometric monitoring of the prototypical symbiotic 
star Z And revealed the highest maximum of its brightness that has ever been recorded by the 
multicolor photometry [U ~8.0, Fig. 1). Around the mid of August the brightness declined to 
f7 ~ 9 and persisted around this level to 2007 January. During this period, rapid photometric 
variation {/S.B ~ Ay ^ 0.06 mag) on the timescale of hours developed (Fig. 3). 

(ii) The SED models and spectral characteristics from the optical maximum can be explained 
by a disk-like structure of the hot active object (Sect. 3.3). The short-term photometric variation 
was produced by the disk (Sect. 3.1). 

(iii) During the optical maximum a mass ejection from the active object was indicated pho- 
tometrically through an increase of the nebular emission (Fig. 2) and spectroscopically by the 
development of complex absorptions on the blue side of the hydrogen and helium lines and a strong 
S"*" jet emission (Fig. 4). 

(iv) High velocity satellite components to Ha and H/3 emission lines developed in the spectrum 
from the optical maximum. Their presence in the spectrum was transient. Our first and last 
detection was on 25/07/2006 and 27/12/2006, respectively. Their spectral properties (Table 2) 
indicated ejection of highly collimated bipolar jets. We summarize their main characteristics as 
follows: 

1. The jets were asymmetrical in their velocities for less than one month from their launching in 
2006 July (i?Vg+/i?Vg- = 1.2 - 1.3). After this episode they became suddenly symmetrical 
until the end of their presence (Fig. 6). 

2. The jets were collimated within an average opening angle of 6°.l for the orbital inclination 
i = 76° (Sect. 3.5). 
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3. The average velocity of jets was 4960 kms ^. The ratio Vjet/'^escape ~ 1 imphes mass of the 
accreting WD of ~ 0.64 M© (Sect. 3.6). 

4. The average outflow rate via jets was Mjet ~ 2 x 10~^(i?jet/l AU)^/^ yr~^, during their 
August-September maximum, which corresponds to the emitting mass in jets, ~ 6 x 
lO-^°(i2jet/lAU)3/2M0 (Sect. 3.8.2). 

5. The total mass released by the jets during their detection was approximated to M^°^^^ 
7.4 X 10-^ M0 (Sect. 3.8.3). 

(v) The short duration of the asymmetry in the jet velocities and evolution in the rapid 
photometric variability resulted from a disruption of the innermost disk. In particular, we pointed 
out a possibility of the radiation-induced warping, which can occur close to the disk's center due 
to the additional source of energy from the thermonuclear burning on the WD surface. For Z And 
parameters, this type of disk instability could occur at all radii R > 0.6 i?© from the accretor 
(Sect. 4.2). 

The authors thank Mario Livio for constructive comments that helped improve the clarity of 
the paper. Dmitri Bisikalo and Drahomir Chochol arc thanked for discussions. Martin Vanko and 
Pavel Schalling are thanked for taking some photometric observations and Miroslav Slechta for 
initial reduction of the spectra taken at the Ondfejov Observatory. This research was supported 
in part by a grant of the Slovak Academy of Sciences No. 2/7010/27 and by the Grant Agency 
of the Czech Republic, Grant No. 205/06/0217 and 205/08/H005. JB acknowledges the support 
from the Pennsylvania State University and the NSF-NATO grant DGE-0312144 and the Marie Curie 
international reintegration grant MIRG-CT-2007-200297. The Hobby-Ebcrly Telescope (HET) is 
a joint project of the University of Texas at Austin, the Pennsylvania State University, Stanford 
University, Ludwig-Maximillians-Universitat Miinchen, and Georg-August-Universitat Gottingen. 
The HET is named in honor of its principal benefactors, William P. Hobby and Robert E. Eberly. 



A. Ionization boundary of jets 

Here we calculate the boundary between the ionized and neutral part of the jet for the case 
that the sources of both the ionizing photons and particles arc located at the central star. The 
boundary is determined by the locus of points, at which the flux of ionizing photons is balanced 
by the rate of ionization/rccombination acts inside the nebula. Thus the equilibrium condition 
between the number of ionizing photons Lph, emitted into the fraction of the sphere AO/47r around 
the direction of the jet, and the number of recombinations, can be expressed as 

Lph— = — / n^npaB{H,T^)4Trr^ dr. (Al) 
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Fig. 9. — Number of photons capable to fully ionize neutral hydrogen and helium as a function of temper- 
ature. The phot on rates were determined by integrating the Planck curve from the ionization limits vq to oo 
(see Eq. (11) in I SkopaTI 1 200 li l . and scaled to the luminosity of lOOOOL©. Solid thin vertical lines limit the 
range of temperatures derived from fluxes of H/3 and He II 4686 during the stable stage of jets (Sect. 3.7). 



It is assumed that the ionizing photons are emitted by the central source spherically symmetrically, 
ael-ffj Je) stands for the total hydrogenic recombination coefficient in the case B, and rip are 
concentrations of electrons and protons and r measures the distance from the central ionizing source 
along the jet. For the fully ionized hydrogen plasma (ng = np) with a mean particle concentration, 
njet, Eq. (Al) can be integrated to give 

^ = ^aB(i^,Te)4t4f (A2) 

In the sense of this relation, the radius of the visible part of the jet, iJjet, represents that of the 
Stromgren sphere. 



B. Number of ionizing photons 

To determine the total number of photons, Lph, capable to ionize neutral hydrogen and helium, 
emitted by the star of temperature T, we integrated the Plank curve from the corresponding 
ionization limits (i.e. 13.598 eV, = 3.29 x 10^^ s"^ and 54.416 eV, = 13.16 x lO^^s-^ for 
hydrogen and helium, respectively) to oo. Figure 9 shows the Lph quantity as a function of the 
temperature for a luminosity of the ionizing source of IOOOOL0. It is worth to note that the 
function Lph{T,H) is slightly decreasing from ~73 000K to higher temperatures, in spite that the 
total energy beyond the ionization limit increases. Therefore, irrespective to a large temperature 
changes between 50 000 and 125 000 K, the number of the hydrogen ionizing photons varies only 
within 15% of the maximum value at 73 000 K at a given luminosity (see Fig. 9). 



-29- 



REFERENCES 

Belloni, T., Klein-Wolt, M., Mendez, M., van der Klis, M., & van Paradijs, J. 2000, A&A, 355, 271 
Bessell, M. S. 1983, PASP, 95, 480 

Bisikalo, D. V., Boyarchuk, A. A., Kilpio, E. Yu., Tomov, N. A., Tomova, M. T. 2006, Astron. Rep., 80, 722 

Brocksopp, C, Sokoloski, J. L., Kaiser, C, Richards, A.M., Muxlow, T.W.B., & Seymour, N. 2004, MNRAS, 
347, 430 

Burmeister, M., & Leedjarv, L. 2007, A&A, 461, L5 

Crocker, M. M., Davis R. J., Eyres S. P. S., Bode M. F., Taylor A. R., Skopal A., Kenny H. T., 2001, 
MNRAS, 326, 781 

Fekel, F. C, Hinkle, K. H., Joyce, R. R., & Skrutskie, M. F. 2000, AJ, 120, 3255 

Fernandez-Castro, T., Gonzalez-Riestra, R., Cassatella, A., Taylor, A. R., & Seaquist E. R. 1995, ApJ, 442, 

366 

Formiggini, L. & Leibowitz, E. M. 1994, A&A, 292, 534 
Galloway, D. K., & Sokoloski, J. L. 2004, ApJL, 613, 61 

Gromadzki, M., Mikolajewski, M. Tomov, T., Bellas- Velidis, I., Dapergolas, A., & Galan, C. 2006, Acta 

Astron., 56, 97 

Gurzadyan, G. A. 1997, The Physics and Dynamics of Planetary Nebulae, (Berlin: Springer) 
Hirth, G. A., Mundt, R., Solf, J., & Ray, T. P. 1994, ApJ, 427, L99 
Horn, J., Kubat, J., Harmanec, P., et al. 1996, A&A, 309, 521 

lijima, T. 1981, in NATO Advanced Study Institute 69, Photometric and Spectroscopic Binary Systems, 517 

Izumiura, H. 1999, in Observational Astrophysics in Asia and its Future, ed. P. S. Chen, (Yunnan Observa- 
tory: Chinese Academy of Sciences) p. 77 

Karovska, M., CariUi, Ch. L., Raymong, J. C, & Mattel, J. A. 2007, ApJ, 661, 1048 

Kellogg, E., Anderson, C, Korreck, K., DePasquale, J., Nichols, J., Sokoloski, J. L., Krauss, M., Pedelty, J. 

2007, ApJ, 664, 1079 

Lamers, H. J. G. L. M., & Cassinelli, L. P. 1999, Introduction to Stellar Winds (Cambridge: Cambridge 
Univ. Press) 

Livio, M. 1997, in ASP Conf. Ser. 121, Accretion Phenomena and Related Outflows, ed. D. T. Wickra- 
mashinge , G. V. Bicknell, & L. Ferrario, (San Francisco, CA: ASP), 845 

Livio, M. 2004, Baltic Astron., 13, 273 

Livio, M., & Pringle, J. E. 1997, ApJ, 486, 835 

Livio, M., Pringle, J. E., & King, A. R. 2003, ApJ, 593, 184 

Lizano, S., Heiles, C, Rodriguez, L. F., Koo, B.-C, Shu, F. H., Hasegawa, T., Hayashi, S., & Mirabel, I. F. 
1988, ApJ, 328, 763 

Lopez-Martin, L., Cabrit, S., & Dougados, C. 2003, A&A, 405, LI 



-30- 



Mikolajewska, J., & Kenyon, S. J. 1992, MNRAS, 256, 177 
Mikolajewska, J., & Kenyon, S. J. 1996, AJ, 112, 1659 
Miirset, U., & Nussbaumer, H. 1994, A&A, 282, 586 
Nagae, T., Oka, K., Matsuda, T., et al. 2004, A&A, 419, 335 
Namouni, F. 2007, ApJ, 659, 1505 
Nussbaumer, H., & Vogel, M. 1987, A&A, 182, 51 

Osterbrock, D. E. 1989, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei (San Francisco, CA: 
Freeman) 

Parimucha, S., & Vanko, M. 2005, CoSka, 35, 35 

Pringle, J. E. 1993, in Astrophysical Jets, ed. D. Burgerela, M. Livio & C. O'Dea (Cambridge: Cambridge 
Univ. Press), 1 

Pringle, J. E. 1996, MNRAS, 281, 357 

Pringle, J. E. 1997, MNRAS, 292, 136 

Ramsey, L. W., et al. 1998, in Proc. SPIE, 3352, 34 

Schmid, H. M., & Schild, H. 1997, A&A, 327, 219 

Schmid, H. M., Krautter, J., Appenzeller, I., et al. 1999, A&A, 348, 950 
Schmid, H. M., Kaufer, A., Camenzind, M., et al. 2001, A&A, 377, 206 
Shahbaz, T., Livio, M., Southwell, K. A., & Charles, P. A. 1997, ApJ, 484, L59 

Shapiro, S. L., & Tcukolsky, S. A. 1983, Black Holes, White Dwarfs, and Neutron Stars (New York: Wiley) 

Skopal, A. 1998, A&A, 338, 599 

Skopal, A. 2001, A&A, 366, 157 

Skopal, A. 2003, A&A, 401, L17 

Skopal, A. 2005, A&A, 440, 995 

Skopal, A. 2006, A&A, 457, 1003 

Skopal, A. 2007, New Astron., 12, 597 

Skopal, A., Bode, M.F., Crocker, M. M., Drechsel, H., Eyres, S. P. S., & Komzik, R. 2002, MNRAS, 335, 
1109 

Skopal, A., Pribulla, T., Vahko, M., Vehc, Z., Semkov, E., Wolf, M., & Jones, A. 2004, CoSka, 34, 45 
Skopal, A., & PribuUa, T. 2006, Astron. Tel., 882 
Skopal, A., & Wolf, M. 2006, Astron. Tel., 930 

Skopal, A., Vittone, A. A., Errico, L., Otsuka, M., Tamura, S., Wolf, M., & Elkin, V. G. 2006, A&A, 453, 
279 

Soker, N. & Regev, O. 2003, A&A, 406, 603 



-31 - 



Soker, N. & Lasota, J.-P. 2004, A&A, 422, 1039 
Sokoloski, J. L., & Bildsten, L. 1999, ApJ, 517, 919 
Sokoloski, J. L., & Kenyon, S. J. 2003, ApJ, 584, 1021 

Sokoloski, J. L., Kenyon, S. J., Brocksopp, C, Kaiser, R. C, & Kellogg, E. M. 2004, RevMexAA, 20, 35 
Sokoloski, J. L., Kenyon, S. J., Espey, B. R., et al. 2006, ApJ, 636, 1002 
Stute, M., & Camenzind, M. 2005, A&A, 432, L17 
Stute, M., & Sahai, R. 2007, ApJ, 665, 698 

Skoda, P. 1996, in ASP Conf. Ser. 101, Astronomical Data Analysis Software and Systems V, ed. G. H. 
Jacoby, & J. Barnes, (San Francisco, CA: ASP), 187 

Taylor, A. R., Seaquist, E. R., & Mattel, J. A. 1986, Nature, 319, 38 

Tomov, T., Munari, U., & Marrese, P. M. 2000, A&A, 354, L25 

Tomov, N. A., Taranova, O. G., & Tomova, M. T. 2003, A&A, 401, 669 

Tomov, N. A., Tomova, M. T., & Bisikalo, D. V. 2007, MNRAS, 376, L16 

Tomov, N. A., Tomova, M. T., & Bisikalo, D. V. 2008, MNRAS, 389, 829 

TuU, R. G. 1998, in Proc. SPIE, 3355, Optical Astronomical Instrumentation, ed. S. D'odorico (Bellingham: 
WA: SPIE), 387 

Wijers, R. A. M. J., & Pringle, J. E. 1999, MNRAS, 308, 207 

Woitas, J., Ray, T. P., Bacciotti, F., Davis, C. J., & Eisloffel, J. 2002, ApJ, 580, 336 



This preprint was prepared with the AAS IM^^X macros v5.2. 



